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Abstract 

The generalized Dirac equation of the third order, describing par- 
ticles with spin 1/2 and three mass states, is analyzed. We obtain the 
first order generalized Dirac equation in the 24-dimensional matrix 
form. The mass and spin projection operators are found which ex- 
tract solutions of the wave equation corresponding to pure spin states 
of particles. The density of the electromagnetic current is obtained, 
and minimal and non-minimal (anomalous) electromagnetic interac- 
tions of fermions are considered by introducing three phenomenologi- 
cal parameters. The Hamiltonian form of the first order equation has 
been obtained. 

PACS: 03.65Pm; ll.lOEf; 12.10Kt 

1 Introduction 

The standard model (SM) of electroweak interactions does not explain the 
mass spectra of leptons and quarks, and contains many free parameters. 
One of the most fundamental problems is the mass generation of fermions. 
Therefore, it is important to investigate the unified description of quarks and 
leptons and understand the origin of fermion masses. The mass generation 
mechanism is beyond the SM and is needed for the deeper understanding of 
lepton-quark families. 

We consider here the toy model based on the generalized Dirac equation 
of the third order (see [l]^ possessing flavor 5*3 symmetry, for particles with 
three nonzero mass states. One can speculate that this equation may be ap- 
plied for unifled description of lepton-quark generations (see [2], [3], [5] 

^In the work [IJ, authors investigated the case of two nonzero masses and one zero mass 
for the unified description of e, /x-leptons and neutrino. 



for the case of two mass states). According to the Barut approach [2], such 
equations may be treated as effective equations for partly "dressed" fermions, 
and such schemes represent the non-perturbative approach to quantum elec- 
trodynamics. 

It should be noted that higher derivative (HD) field theories have been 
investigated already. Some examples include generalized electrodjTiamics [S] , 
renormalizable gravity theory in four dimensions [7] and others. Therefore, 
it is motivated to investigate HD field theories in four dimensions. 

In this paper, we represent the third order (in derivatives) equation in 
the form of the first order generalized Dirac equation, obtain solutions in the 
form of projection operators and consider minimal and non-minimal electro- 
magnetic interactions. 

The paper is organized as follows. In Sec. 2, we define the generalized 
Dirac equation of the third order for fermions, obtain projection operators 
extracting solutions of the HD equation, and analyze a propagator of fields. 
The first order generalized 24-component Dirac equation is derived in Sec. 
3. It is proven that the mass spectrum of the 24-component matrix equation 
obtained and the third order 4-component equation is the same. We obtain 
mass and spin projection operators in Sec. 4. These projection operators ex- 
tract pure spin states of particles. Sec. 5 is devoted to introducing minimal 
and non-minimal electromagnetic interactions of fermions. In this section, we 
find also the electromagnetic current density. We postulate the matrix equa- 
tion with three phenomenological parameters characterizing anomalous elec- 
tromagnetic fermion interactions and derive the quantum-mechanical Hamil- 
tonian. A conclusion is made in Sec. 6. In appendixes, we obtain useful 
relations of 24 x 24-matrices of the wave equation and spin matrices. 

The system of units h = c = 1 is chosen and euclidian metric is used (see 

El). 

2 Field Equation of Third Order for Fermions 

We postulate the third order (in derivatives) field equation describing fermions: 

ilfid^ + mi) {-f^du + m2) {jada + m^) ^(x) = 0, (1) 

where = d/dx^ = {d/dxm,d/d{it)), ip{x) is a bispinor; mi, m2, are 
masses of fermions. As usual, repeated indices imply a summation. The 
Dirac Hermitian matrices 7^ obey the commutation relations 7^7jy + 'Ju'j^ = 



2 



25^1,. It is obvious, that there are three independent solutions of Eq. (1) 

corresponding to masses m.i {i = 1,2,3). The bispinor ip^x) ~ ip~^{x)^4^ ("0"*" 
is the Hermitian-conjugated bispinor) satisfies the equation 

7jj{x) (jiu*dn - mi) (ov^^ - 7712) (ja*da - ms) = 0, (2) 

where the derivatives d ^ act on the left-standing function. The discrete 
symmetries of the wave function ip^x) directly follow from the Dirac theory. 
Eq. (1) for positive energy, in momentum space, is given by 

(ip + mi) (ip + 1112) (ip + 1713) ijj(p) ^ 0, (3) 

where p — j^Pfj,, p^ = (p, ipo)- For antiparticles one should make the re- 
placement p — > —p. Prom Eq. (3), one obtains the dispersion equation as 
follows: 

(y + ml) (y + ml) + ml) = 0, (4) 

The projection operator extracting solutions to Eq. (3), corresponding to 
positive energy, is given by 

^ ^ (-ip + mi) {-ip + 7712) {-ip + ma) 
2[mim2m3 - p'^{mi + m2 + m2)] 

It is easy to verify that the projection operator (5) obeys the necessary equa- 
tion = A+. It should be noted that the projection operator A_|_ is the 
solution to Eq. (3) for any momentum p obeying the dispersion relation 
(4), i.e. the momentum p is not specified here. The projection operator A_|_ 
corresponds to positive energy, and at the replacement p — > —p, to negative 
energy of particles. 

In the momentum space the propagator of the HD field is 

^(P) = 77 



{ip + mi) {ip + m2) {ip + ms) 

(6) 

ip + mi){—ip + m2){—ip + m^) 



{p^ + mf) {p^ + ml) {p^ + ml) 
The propagator (6) can be represented as a sum of Dirac propagators 

^ A{-ip + mi) _ B{-ip + m2) C{-ip + mg) 
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where the coefficients A, B, and C are given by 



A 



1 



B 



1 



(m2 — mi)(m3 — mi) ' 



{■1712 - "T-Ol^is - "^2) ' 



(8) 



C 



1 



(777,3 - m2){m^ - mi) 



Eq. (7) is vahd only in the case 777i 7^ 7712 7^ 7713, i.e. masses of the field 
are not degenerated. In the case 7773 > 7772 > 777i, the coefficients A, B, and 
C are positive, and the second propagator in Eq. (7), corresponding to the 
mass 7772, has the "wrong" sign (— ). This means that the HD field possesses 
two physical states with masses 7771, 7773, and one Weyl ghost with the mass 
7772. The ghost gives the negative contribution to the energy [3], [9], and, as a 
result, the Hamiltonian is not bounded from below. Therefore, to formulate 
the quantum field theory one needs to introduce the indefinite metrics [3]. 
As the propagator (7) is a sum of Dirac propagators, the renormalization 
of the theory can be performed with the help of the standard procedure of 



3 24-Component First Order Relativistic Wave 
Equation 

For a convenience, we rewrite Eq. (1) in the form 



where = d'^ — dl (777 = 1, 2, 3, = d/dt). The coefficients 

a = 777i 7772 + 777 i7773 + 77727773, 6 = 771i + 7772 + 7773, C = 777.177727773 (10) 

are invariant under the permutation symmetry 5*3. 

Now we reformulate the third order Eq. (9) in the form of the first order 
relativistic wave equation. It is useful for Lagrangian formulation of the 
theory, obtaining the conserved currents, and different calculations. Let us 
introduce the 24-dimensional function 



QED. 






(11) 
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where A = 0,0, /i, and 

^'o(x) = ^{x), ^o(x) = ^{x) = -dltPix), 

a ^ 

(12) 



The wave function represents, in Eq. (12), the direct sum of bispinor 
ip{x), bispinor il){x), vector-bispinor %l)^[x), and transforms under the Lorentz 
group as 

[(1/2, 0) © (0, 1/2)] © [(1/2, 0) © (0, 1/2)] © {(1/2, 1/2) © [(1/2, 0) © (0, 1/2)]}. 

To obtain the first order relativistic wave equation from Eq. (9), we introduce 
the elements of the entire algebra e^'^ (see, for example, [10]) obeying the 
multiplication rule and having matrix elements as follows: 

£ ' £ ' = OabS , [£ )^^=OmaOnB, (lo) 

A,B,M,N = 0,0,/i (/i = 1,2,3,4). The matrix e^^'^ consists of zeros and 
only one element is unity, where row M and column cross. 

Using Eq. (11)-(13), Eq. (9) takes the form of the first order equation 

\ J AB 

(14) 

+ f + + fe^^'^) vl/^(x) = 0. 
V a J AB 

It is implied that 7-matrices act on the bispinor indexes. There is a sum- 
mation over repeated indices in subscripts and superscripts of the matrices. 
Now we introduce the matrices 

= l^e^'^ + 5°'^ - eO-'^) © /4 + + © 7m, (15) 

M = -Po + bPo + ^P^, (16) 
a 

Po=£0'0®/4, p^ = e°'°©/4. Pi =£'^'^©14. (17) 

Unit four-dimensional matrix I4 acts on bispinor indexes and matrices e^'^ 
{M,N = 0,0, fi) act on scalar and vector indexes. It should be mentioned 
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that matrices Pq, Pq, Pi are projection matrices obeying the relation P^ = P. 
The M is the diagonal and Hermitian matrix, M = M"*", but the matrix 
is not Hermitian matrix. 

Taking into account Eq. (11)-(13), Eq. (14) becomes 

{T^d^ + M)^{x) = 0. (18) 

Eq. (18) obtained, is the relativistic 24-component wave equation of the first 
order and is convenient for different applications. 

Now, we derive the mass spectrum of Eq. (18). In the frame of references 
where a particle is at the rest, Eq. (18), in the momentum space, becomes 

(por4-M)^(po) = 0, (19) 

where Pq is the energy of a particle. Here we use one particle (quantum me- 
chanical) interpretation of Eq. (18) and do not discuss the second quantized 
theory. Masses of particles (see [12]) are given by the relation = A~\ 
where Aj are eigenvalues of the matrix M'^T^ or rrii = Aj, where Aj are 
eigenvalues of the matrix MFj^. From Eq. (15), (16), one may find expres- 
sions 

= -Po + jPo + -Pi, (20) 
c a 

K ^ M-IF4 = { -e''' + le'^' - -e'A 0h+( + -e'A ® 74, (21) 
\a b c y Vc c / 

With the help of Eq. (13), it is easy to verify that the matrix K obeys 
the equation as follows: 

cK^Ie ® 74 - aK^ + hKh ® 74 - A = 0, (22) 

where 

/e ® 74 = + + ^^''') ® 74, (23) 

A = Po + + ^''^ ® h. (24) 

The matrix A is the projection matrix, so that A^ = A. It follows from 
Eq. (22), that as the matrix /g ® 74 possesses eigenvalues ±1, the matrix 
A possesses eigenvalues 1, eigenvalues of the matrix A = {Aj}, obey the 
equation 

± cA^ - aA^ ± 6A - 1 = 0. (25) 
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As a result, taking into consideration the relation rrii — the mass spec- 
trum of Eq. (18) is given by solutions of the equation 

T bml + ami T c = 0. (26) 

Positive solutions to Eq. (26), rrii, are connected with the coefficients a, 
6,and c, possessing the flavor 5*3 symmetry, by relationships (10). So, the 
mass spectrum of Eq. (18) and Eq. (9) is the same. 

As the diagonal mass matrix M is not proportional to unit matrix, it 
is convenient for some applications to get the equivalent form of Eq. (18). 
Multiplying Eq. (18) by the non-singular matrix 



—r^o + - 
ab a 

and introducing the new matrix 



Eq. (18) transforms to the standard form of the relativistic wave equation 

(T,d, + m)<i/{x)=Q. (28) 
We took into consideration the equalities 

N M = -/24, m = 



a a 171x1712 + 171x1713 + 17121713 

where I24 is the unit 24 x 24-matrix, and m is a reduced mass. Notice that 
the wave function ^{x) remains the same (Eq. (11)) in Eq. (18), (28). 



4 Mass and Spin Projection Operators 

In the momentum space Eq. (28) becomes 

{ip + m) ^(p) = 0, 



(29) 



where p — p^T. One may verify, with the help of Eq. (A6)-(A7) (see Ap- 
pendix A), that the p obeys the matrix equation as follows: 



p{p + m ] c p — m p [pa + cj—mp 



0, 



(30) 
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where d = m\m\ + m\m\ + m\m\. From Eq. (30), we find the projection 
operator 

„ pip + im) [c^p'^ — m?p^ {p^d + c^) — m^p^] . . 

+ " 2m6 (p6 _p2^_2c2) • ^ ^ 

It is not difficult to verify, with the help of Eq. (30) that the the projection 
operator (31) obeys the equation 11^ = n+. To find the projection operator 
corresponding to negative energy, n_, one has to make the replacement p — *• 
—p in Eq. (31). Every column of the matrix 114. is a solution to Eq. (29), 
and after acting on 24-dimensional vector x, one finds the solution to Eq. 
(29), \l/(p) = n+x. The matrix (31) also represents the density matrix for 
the impure spin state. 

To obtain the spin projection operators, we consider the generators of the 
Lorentz group in the 24-dimensional representation which are given by 

J;.. = I I ®h + h®J^^J^\ (32) 




where [TH] 



= e^,- _ e-./^, (33) 

J^J'^ = \ {1,1. - 7.7.) • (34) 

Unit 6-dimensional matrix Jg is defined in the space (0, 0, /x), and the genera- 
tors of the Lorentz group in four- dimensional vector and bispinor spaces are 
J^J ^ J^J'^\ respectively. It is easy to verify that generators (32)-(34) satisfy 
the commutation relations 

In our euclidian metric, the antisymmetric parameters of the Lorentz group 
^mn {m,n = 1,2,3) are real, and cOmi are imaginary. The relativistic form- 
invariance of Eq. (28) follows from the commutation relation [T2] 



hfJiy - hiy^/i, (36) 



which is valid for matrices (27), (32). 
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The operator of the spin projection on the direction of the momentum p 
is given by (see, for instance, [TU]) 



I 

= ~ly~\^a.bcVaJbc-, (37) 



where |p| = ^Jp'l + + pI, €abc {a,b,c = 1,2,3) is an antisymmetric tensor 
Leyy-Civita. With the help of Eq. (36), one may verify that the operator 
of Eq. (28) in the momentum space (^F^p^ + m) commutes with the spin 
operator (37): [iT^p^ + m,ap] = 0. This guaranties the existence of the 
common solutions of the Eq. (29) in the momentum space, "^^^^Ip), and the 
equation 

where Sp = ±1/2. With the aid of Eq. (13), (B5) (see Appendix B), one can 
verify that the matrix equation 

2 ^\ f 2 9 



is valid. This equation allows us to construct two projection operators cor- 
responding to spin 1/2 and 3/2. Indeed, the representation of the Lorentz 
group corresponding to the 24-dimensional wave function (11), contains also 
the spin-3/2 component. It is important to show that the wave equation (28) 
(or (18)) describes only pure spin-1/2 states. For this purpose, we consider 
both projection operators corresponding to spin 1/2 and 3/2. Exploring the 
method [H], one obtains the projection operators 

1 Z . n / , 9\ „ _1 / , 3\ / , 1 



P.U^ = ± 2 j [< - i j . ^«/. = [a, ± 2 j ["i - -J (40) 

extracting spin projections ±1/2 and ±3/2. One can verify that = 
{sp = ±1/2, ±3/2). With the help of Eq. (B6), we find the relations" 

P±3/2P = pP±3/2 = 0, P±3/2n± = n±P±3/2 = 0, (41) 

but P±i/2n± 7^ 0. Therefore, the wave equation (28) (or (18)) describes only 
pure spin-1/2 states. Eq. (39), (40) lead to the equation 

apP,^ = SpPs^. (42) 



9 



Eq. (42) shows that the wave function "^spip) — PspX- (x being an arbitrary 
nonzero 24-dimensional column) is the solution of Eq. (38). The projection 
operator corresponding to common solutions to Eq. (29), (38) is given by 

P±i/2(p) = n+p±i/2. (43) 

The projection operator (43) extracts states with positive energy of particles 
and spin projections ±1/2. The p±i/2 is also the density matrix for pure spin 
states (sp = ±1/2). 



5 Non-Minimal Electromagnetic Interactions 
of Fermions 

Now, we obtain the density of the electromagnetic current. For this purpose, 
Eq. (2) is rewritten as follows: 



(44) 



Following the well-known procedure, multiplying Eq. (9) (from left) by ipix), 
and Eq. (44) (from right) by ip{x), and adding them, one finds 



ipix)j^,di,dlip{x) + (d^dlip{x)) j^ip{x) 



+ a 
+b 



(45) 



\p{x)dlip{x) - (aXx)) ip{x) 



0. 



Identifying the left side of Eq. (45) with (— ia)9juj^(x), we obtain the electric 
current density 



6 

ju,{x) = i^}){x)'^^il){x) +i- ^l){x)d^^l){x) - {d^i){x)) i){x) 



+ -[^{.x)-i^dlij{x) + (dlil){x)\ 'y^.'^ix) - (d^ipix)) 7^a^^(a;) 

(X 



(46) 



so that the conservation of the four- vector current density df^jfj_{x) = holds. 
The current jfj,{x) consists of the usual Dirac current j^{x) — i'il){x)^^'^{x) 
and additional convective terms. 
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Substituting the derivatives in Eq. (18), = — ieA^ {A^ is 

the four- vector potential of the electromagnetic field), one can obtain the 
minimal interaction with an electromagnetic field. Now we introduce the 
non-minimal electromagnetic interaction by considering the matrix equation 



(47) 



where jF^j, = d^A^, — d,^A^ is the strength of the electromagnetic field, and 
kq, kq, Ki are parameters which characterize anomalous electromagnetic in- 
teractions of fermions, and 



r = r r 



r r 



(48) 



The projection operators Pq, Pq, Pi obey the relations: PqPq — PqPi — 
PqPi = 0, Pq + Pq + Pi = /24. Eq. (47) is the relativistic wave equation 
for interacting fermions which is form-invariant under the Lorentz transfor- 
mations. The tensor form of Eq. (47) follows from Eq. (11), (15)- (17), 



(49) 



Df,ijj{x) + iJ^^^ulv Koip{x) + Koiplx) 



-6^y - iKiJ^^u ) i)u{x) = 0, (50) 



D^xlj^{x)+hjj{x)^0. (51) 

Eq. (49)-(51) can be used for the phenomenological apphcations in the case 
of anomalous electromagnetic interactions of fermions. The system of Eq. 
(49)-(51) defines bispinors ip{x), ip{x), and vector-bispinor ipu{x)- Eq. (50) 
may be considered as a matrix equation for vector-bispinor ipiy{x). Finding 
ip{x), and vector-bispinor ipu{x) from Eq. (51), (50) and substituting them 
into Eq. (49), one may obtain an equation for bispinor ip{x). This equation 
contains the interaction with the anomalous magnetic moments of particles. 
To figure out the physical interpretation of introduced constants Kq, Ki, Kq, 
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one needs to consider the non-relativistic limit of the theory. We leave this 
for further investigations. 

Let us consider the Hamiltonian form of Eq. (47). This form is convenient 
for determination of dynamical variables of the wave function "^(x) and for 
second quantization of fields under consideration. Eq. (47) can be rewritten 
as 



2 

(52) 

With the aid of Eq. (13), one may verify that the matrix r4 obeys the 
equation as follows: 

d = A, (53) 

where the projector operator A is given by Eq. (24). To separate the canon- 
ical and non-canonical parts of the wave function, we define the projection 
operator 

n = 1 - A, (54) 

where 1 = I24 is the unit 24 x 24-matrix. Projection operators A, 11 obey 
equations 

n + A = 1, nA = An = 0, a^ = a, = n. (55) 

Introducing the auxiliary wave functions 

0(x)=A*(x), x(^) = n*(x), (56) 

so that (j){x) + x{^) = ^(2^)) multiplying Eq. (52) by r|, and taking into 
account Eq. (53), one obtains 

idt(l>(x) = e A0(x) + rl (r^L>^ + (7 + m) [0(x) +x{x)], (57) 

where we have introduced the operator 

cr = ^y^tMu^tMu [kqPo + /^o-Po + «i-Pi) • (58) 
Acting on Eq. (52) with the operator 11, and using the equations 

nr4 = 0, UM = yU, nr^n = o, 
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we obtain the equation for auxiliary wave function 

^X{x) + nTmDm(j){x) + Ua Mx) + xi^)] = 0. (59) 

Expressing the auxihary wave function xi^) from Eq. (59), and replacing it 
into Eq. (57), one obtains the equation in the Hamiltonian form: 

idt(f){x) — 'H(f){x), 

H^eAo + rllrmDm + a + M^ l-(^^ + Ua^ H (r^D^ + a) 

The wave function (j){x) entering Eq. (60), possesses 12 components, and 
the auxiliary wave function x(x) also has 12 components. Indeed, for the 
relativistic description of spin-1/2 fields (four components for each field) with 
three mass states, 12 components are needed, and, therefore, the dynamical 
wave function (j){x), possesses 12 components. 

At the particular case, for point-like particles, without anomalous inter- 
actions, (7 = 0, one finds the Hamiltonian from Eq. (60) as follows: 



(61) 



6 Conclusion 

We have considered the generalized Dirac equation of the third order, describ- 
ing fermions possessing the permutation symmetry 5*3. The 24-dimensional 
matrix form of the first order was obtained which is convenient for differ- 
ent applications. Although the dimension of the matrices is high (24 x 24), 
we have found minimal equations for matrices. This allowed us to obtain 
the projection operators extracting states with definite energy and spin pro- 
jections. The density of the electromagnetic current was obtained, and we 
have introduced three phenomenological parameters characterizing anoma- 
lous non-minimal electromagnetic interactions of fermions. It is of interest to 
clear up the physical interpretation of constants introduced. This and other 
questions we leave for further investigations. 

Dynamical and non-dynamical variables were separated and quantum- 
mechanical Hamiltonian was derived. Fist order matrix wave equations ob- 
tained are convenient for second quantization and for calculations of different 
electro dynamic processes. 
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Non-trivial differences between this three mass case and the previously 
studied two mass generalized Dirac equation p] is in the dimension of matri- 
ces of wave equation and corresponding minimal equations. We mention that 
in this three mass case it is not trivial to obtain the Hermitianizing matrix 
r] in our 24- representation case and corresponding Lagrangian. 

The HD field theory considered includes the ghost state which requires to 
introduce indefinite metrics. It should be noted that this negative attribute 
appears in many HD field theories. 

One of important unanswered questions is the mechanism of mass gener- 
ations of fermions (see [13] and references there). Here, we have introduced 
three independent parameters which specify the masses and do not explain 
their appearance. It seems interesting to follow the Barut [2] idea about 
the electromagnetic nature of fermion masses. This, however, requires the 
further investigations of the approach considered. 



For a convenience, we represent the matrix p of Eq. (29) as follows 



Appendix A 



p = /(o) ®P + /(I) ® h 



where 





a'' ab 
Using Eq. (13) one may prove simple relationships: 

1(0)1(1)1(0) = 0, /fi) = J(i)/(o)/( 




{A3) 




1(2) = P,,Pye' 



'2 2t 
(2) =P H2), 






1(2)^(1) = -a/ ^(2), ^(o)^(i)ho) = «(/^ - l)P^^(o), 

^a)Ao)^(i) + 1(^)^(0)^(1) = - 
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One may check, with the help of Eq. (13),(A4), that the minimal equation 

/ + {2a - + a{a - 2(5)p^p^ - a^pYp = 0. (^5) 

is valid. It should be noted that the dispersion equation (4) can be rewritten 
as 

- (2a - l)my + a{a - 2/?)mV + a^/^V = 0. {A6) 
With the help of Eq. (A6), Eq. (A5) may be represented as follows: 

Although the p is 24 x 24-matrix it obeys the simple minimal equation. 

Appendix B 
The spin operator (37) can be written as 

(7p = (7«(8)74 + /6(H)<7i'/'\ (Bl) 

where 



c,. 



p — 2|p| "''^P" ''^ ' P ~ 2|p| "''^P" '"^ ■ ^ ^ 



The operators a^^\ CTp^'^^ correspond to spin-1 and spin-1/2 representations 
of the Lorentz group. Using Eq. (33), (34), we obtain 

so, that the spin operators obey the simple matrix equations 

(.<■) - 1) (.<■) + 1) = 0. - i) (.<V^. + 1) = 0. 

From Eq. (Bl), (B3), one finds 

= ^ + K'0'®/4 + 2aW®a(V2). (55) 
Using Eq. (B3), (B5), we obtain the useful relation 

^2 1 



p = 0. (56) 
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